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Abstract 

We present an analytic study of the physics of the glasma which is a strong classical 
gluon field created at early stage of high-energy heavy-ion collisions. Our analysis is 
based on the picture that the glasma just after the collision is made of color electric 
and magnetic flux tubes extending in the longitudinal direction with their diameters 
of the order of 1/Q S (Qs is the saturation scale of the colliding nuclei). We find that 
both the electric and magnetic flux tubes expand outwards and the field strength 
inside the flux tube decays rapidly in time. Next we investigate whether there exist 
instabilities against small rapidity-dependent perturbations for a fixed color config- 
uration. We find that the magnetic background field exhibits an instability induced 
by the fluctuations in the lowest Landau level, and it grows exponentially in the 
time scale of 1 /Q s . For the electric background field we find no apparent instability 
while the possible relation to the Schwinger mechanism for particle pair creations 
is suggested. 



1 Introduction 



There is a missing link in our current understanding of the event evolution in 
ultra-relativistic heavy-ion collisions. Before the collision, the projectile nuclei 
at very high energies are described in the framework of the Color Glass Con- 
densate (CGC) [1]. Once the quark-gluon plasma (QGP) is formed in local 
thermal equilibrium, its time evolution may be described with hydrodynamics 
[2]. Obviously here it is a big theoretical challenge how the local equilibrium 
state is established from the non-equilibrium initial condition given by the 
CGC. In particular, ideal hydrodynamic simulations seem to require that, at 
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RHIC, thermalization should be realized within a very short period r < 1 
fm/c, which was a surprise from the viewpoint of the standard theoretical 
estimates. 

Theoretical efforts involve the detailed analyses of the plasma instability (in 
particular, the Weibel instability) [3] which takes place when the hard par- 
ticles having very anisotropic distribution in the momentum space interact 
with soft gauge fields (See Ref. [4] for other scenarios and more references). 
Although many new interesting phenomena have been found in these stud- 
ies, they are based on the kinetic equation applicable for r <: l/Q s when 
individual particles having transverse momenta greater than Q s are formed 
out of fields. The plasma instability scenario, therefore, should be carefully 
examined at earlier stages r ^ 1/Q S , where the system is so dense that the 
field description may be more appropriate. Yet, it is claimed that there is 
a certain intermediate regime where both the Boltzmann and classical field 
equations give almost the same description for the system [5]. It is therefore 
very likely that the plasma instability scenario itself may be extended to such 
earlier times r ^ l/Q s where the kinetic theory should be turned over to the 
classical Yang-Mills equation; thus we shall investigate the instability problem 
in the classical Yang-Mills description. 

The dense pre-equilibrium system appearing between the first impact and the 
equilibrated QGP is recently named Glasma [6]. The glasma, produced from 
the CGC initial condition, is still well-described by strong coherent Yang- 
Mills field, which can be treated as a weak-coupling system. However, unique 
to the glasma is that it has strong time dependence which is believed to lead 
eventually to the QGP. 

If there are unstable modes in the glasma, it will help the matter evolve to an 
isotropic/thermalized state more rapidly. In fact, instability in the glasma is 
already found in the lattice numerical simulation [7]. Properties of the unstable 
glasma turned out to be similar to those of the Weibel instability in expanding 
geometry [8] (see also Ref. [9] for more recent results). On the other hand, the 
instability in the glasma is analytically less understood although there are 
several works about the initial fluctuation of the glasma [10] and its possible 
instability [11,12,13,14,15]. In this paper, we shall present a detailed analytic 
study of the dynamics and instability of the glasma. 

Our theoretical framework is based on the following picture. First of all, we 
look for instabilities of the glasma in a single event to which one fixed color 
configuration is assigned from the initial distribution of the CGC. Namely we 
presume that thermal equilibrium should be achieved for each event without 
taking the average over the distribution of the initial condition. The typical 
initial configuration specified by the CGC is dense color charges randomly 
populated on the two-dimensional transverse plane with its coherence length 
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given by 1/Q S - After the collision, this configuration induces strong electric 
and magnetic fields in the longitudinal direction gE z ~ gB z ~ Q 2 S) while the 
transverse components are vanishing E l = B l = 0. Therefore, the glasma just 
after the collision can be characterized by flux tubes of strong color fields 
extending in the longitudinal direction with their transverse size being 1/Q S . 
This configuration is independent of the rapidity coordinate, provided that the 
glasma is expanding longitudinally at the speed of light. We will investigate the 
time evolution of this flux-tube configuration. Notice that this boost-invariant 
glasma cannot thermalize because it does not allow longitudinal momentum. 
It is necessary to introduce the rapidity dependence to the glasma, which is 
only possible through rapidity-dependent fluctuations. We will analytically 
examine the stability of the rapidity-dependent fluctuations, and find indeed 
that there exists an instability in the magnetic background! 1 1 

The present paper is organized as follows: In the next section, we will define 
more precisely our theoretical framework, and discuss the properties of boost- 
invariant glasma as a collection of longitudinal flux tubes. Then, in Sect. 3, 
we will show a detailed analysis of the rapidity-dependent fluctuations in the 
boost-invariant background. Summary and discussion will follow in the last 
section. 



2 Dynamics of the boost-invariant glasma 



At high energies, heavy-ion collisions can be viewed as the CGC-CGC colli- 
sions [17]. Namely, the incident nuclei, right-moving 1 and left-moving 2, are 
highly Lorentz contracted and the valence and large- x degrees of freedom of 
each nucleus can be represented as static color source, px^{x±)i which has 
random distribution on the transverse plane. The small- a; gluons are treated 
as the radiation field generated by these color charges. Thus, the collision is 
described by the Yang- Mills equation [D^, F^ v ] = J u with the color current 
J v = 5 u+ 5(x-)p 1 (x±) + 5 u -5(x + )p 2 (x ± ), where x ± = (x° ± x 3 )/V2. Notice 
that the current J M is nonzero only on the light-cone axes x = 0. Thus, the 
dynamics of the glasma in the forward light cone x > (after the collision) 
is simply described by the source free Yang-Mills equation [D^, F^ u ] = with 
appropriate boundary conditions on the light-cone axes. It is only through 
this boundary conditions where the information of the colliding CGC's enters. 
In fact, as we will discuss later, the initial condition supplied by the CGC is 
quite unique, and greatly affects the dynamics of the glasma. 

In order to describe the glasma expanding in the longitudinal direction at 
the speed of light, it is convenient to introduce the T-r] coordinates defined 

1 Similar conclusion is recently obtained in Ref. [16] independently of ours. 



3 



by r = \ / 2x + x~ and f] — \ ln(x + /a; ) . In these coordinates, the source free 
Yang-Mills equation [D^^F^] = is rewritten separately for v = r,rj,i as 



r[D T ,^F TV ]-[D t ,F lv ]=0, (1) 
[D v ,-F T ^-[D l ,rF tT } = 0, (2) 

T 

- [.D T , rF lr ] - -i [£>„, F ir? ] + [£>,-, i^] = . (3) 

T T 

The second equation is the Gauss law constraint, which is obvious if one 
rewrites it as [D v , P v ] + [Di, P l ] = by using canonical momenta defined by 
P r > = 5(TC)/5(d T A v ) = \F TTi and P l = 5(tC) / '5(d T Ai) = -rF ir . When we 
solve the Yang-Mills equation, we have to check if the solution indeed satisfies 
the Gauss law. 

For later convenience, we show here the definitions of electric and magnetic 
field strength k = 1, 2, 3): 

E l = F t0 , B l = -\e ijk Fi k . (4) 



Each component is explicitly represented with respect to quantities in the T-rj 
coordinates as (i — 1, 2) 



-Fj„sinh?7, (5) 

T 

(6) 

1 „ 



E l = 


—F iT cosh f] + 




1 


E z = 


-F 




T 


B l = 


€ij ^Fj T sinh?7 


B z = 


—F±2 . 



T 



(7) 
(8) 



In the rest of this section, we discuss the dynamics of the boost-invariant 
glasma paying attention to the characteristic properties of the CGC initial 
conditions. All the calculation is done in the Fock-Schwinger gauge A T = 
\{x + A~ +x-A+) =0. 



2.1 Boost-invariant glasma 



The glasma created in high-energy collisions expands in the longitudinal di- 
rection at the speed of light. In an idealized situation, the gauge fields after 
the collisions are independent of the rapidity coordinate i] [17]: 
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A v = A v = -T 2 a(r,x ± ) , (9) 
A i = A i = a i (T,x ± ). (10) 

We use calligraphic letters to represent boost invariant glasma which corre- 
sponds to the classical background fields for instability problem discussed in 
the next section. These fields satisfy the following equations 



-d T {T*d T a) -\V h [D^a]] = 0, (11) 
-[T>i, d T ai] - igr[a, d T a] = , (12) 

T 

-d T (Td T ai) - igr 2 [a, [D h a}} - [Dj, J 7 ^] = , (13) 
r 

where the covariant derivative T>i and the field strength Tji are, respectively, 
given by [Pj, a] = dia — ig[ai, a] and Tji = djOti — diaj — ig[otj, ««] . 

These equations are to be solved with the following initial condition 



a(r = 0,x ± ) = - l -j[a\(x ± ),ai(x ± )}, (14) 
a\r — 0, x±) — a\(x±) + a l 2 (x±) , (15) 

where ct\(x±) and a l 2 (x±) are the classical gauge fields of colliding nuclei, 
1 and 2, before the collision. Those are associated with the random color 
sources on x ± = 0, and therefore the initial fields a(r = 0,x±) and c^(r = 
0, x±), too, may be considered as random variables. In the original McLerran- 
Venugopalan model, the color charge distribution on the transverse plane is 
Gaussian random and there is no correlation between the charges at different 
transverse points. However, in reality, the charge distribution is not completely 
random. There is some structure in it. For example, it must be color neutral 
if it is averaged over the transverse area comparable with the nucleon radius. 
In fact, one can show that color neutrality is ensured for the transverse area 
whose extent is larger than 1/Q S [18]. In other words, this implies that the 
classical gauge field is randomly distributed on the transverse plane, but its 
coherence length is at most of the order of l/Q s - 

Just after the collision, only the longitudinal components of the electric and 
magnetic fields are nonzero [20,6] 



E z \ T=0+ = -tg[a\,al], (16) 
B z \ T=0+ =ige ij [a\,a 2 : }. (17) 

Note that these fields are rapidity independent as the gauge fields a\ j2 ( x ±) 
of the nuclei are assumed to be rapidity independent. On the other hand, 
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Fig. 1. The color electric and magnetic flux tubes just after the collision. The trans- 
verse size of the flux tube is of the order of 1/Q S - 

as mentioned above, the coherence length in the transverse plane is ~ 1/Q S - 
Therefore, one finds that the initial configuration of the glasma gauge fields 
is represented as collection of color electric and magnetic flux tubes in the 
longitudinal direction with their diameters of the order of 1/Q S as shown in 
Fig. 1. According to the numerical results by Lappi and McLerran [6], these 
longitudinal fields decay rapidly while the other components are generated as 
r proceeds. 

In the following, for the sake of clarity, we separately discuss two cases with 
(i) a = and «j ^ 0, and with (ii) a ^ and «j = 0, which correspond to 
the magnetic and electric flux tubes, respectively. In fact, the results are valid 
even when both a and «j are nonzero, provided that they have the same color. 

2.2 Magnetic flux tube 

Let us first consider the dynamics of a magnetic flux tube, which is a unique 
object in the glasma. We assume that the magnetic flux tube is isolated and 
there is no longitudinal electric field around it. This situation is simply real- 
ized, from Eqs. (16) and (17), when a\ and a\ have the same color but a\ 
and a J 2 (i ^ j) have different colors from each other, which corresponds to the 
initial condition a(r = 0) = and a*(r = 0) ^ in Eqs. (14) and (15). For 
r > 0, we preserve this condition 

a = 0, ai^O. (18) 

Then we have the field strengths: 

E z = , B z = -F l2 , E i = -F iT cosh rj , B i = e ij F jT sinh T]. (19) 
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Note that the transverse components E l and B l are transformed among them 
under the longitudinal boost. We consider the fields only in the local rest-frame 
(?7 = 0), and then E { (t] = 0) = d T a { , B% = 0) = 0. 

Notice that, by using residual gauge freedom, one can always set [a^a^] = 
so that the magnetic field B z = — iq 2 = —(d x a y — d y ot x ) has the same 
color structure as a*. Let us further assume, as a special case, that and 
B z are directed to a certain color. In other words, we assume here that the 
initial conditions (14) and (15) allow the gauge field configurations which can 
be transformed to the Abelian form with a constant color direct ioEpQ. This 
restricts the dynamics to the Abelian subgroup, and one can ignore all the 
commutators in the equation for c^. 

Since a = 0, the first Yang- Mills equation (11) is trivial, and the other two 
equations (12) and (13), respectively, reduce to 



(20) 

= (2D 

where we have already introduced the Fourier mode <5j(r, k±) with respect 
to transverse coordinates. Note that the factor (<5jj — kikj/k'jj is a projection 
operator to the direction perpendicular to ki. Thus, it is convenient to decom- 
pose a two-dimensional vector <Sj into two parts (parallel and perpendicular 
to ki) 

a t = kiC\\ + e ij k j C ± . (22) 



d T kiOti = 



d T ati H — d T o.i + k ± [5^ — 



Substituting this into Eq. (20), one finds that C\\ is independent of r. Indeed, 
this first term is the remaining gauge freedom and unphysical. We set Cii = 0. 
Then the Gauss law (20) is automatically satisfied and Eq. (21) becomes a 
simpler equation for C±: 

8 2 T C ± + -d T C ± + k 2 ± C ± = , (23) 



2 Incidentally, it is known that the covariantly constant field solution to the source- 
free Yang-Mills equation is generally written in a color factorized form: F^ u = 
(F^ u )n a t a and = (F^ u )x u n a t a where (F^ u ) is independent of x M and n a is a 
constant color vector [19]. This situation should be approximately realized inside of 
the flux tube. 
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which is solved by the Bessel function Jo(\k±\r) . Therefore, the solution to 
Eq. (13) with a = which satisfies the initial condition a- mt (x_|_) is given by 

<*i(T,x ± ) = J ^y 2 e^J (\k ± \r) c$*(fcjL) , (24) 

where the Fourier transform of the initial condition must have the following 
form 

d^\k ± ) = J d 2 x ± e~ ik ^ o^\x x ) = <,jk,fik.) . (25) 

Notice that the series expansion of Jo{\k±\r) has only terms with even powers 
of r. This is consistent with the result of r-expansion developed in Ref . [20] . It 
is also important to notice that the Bessel function Jq(z) is an oscillating func- 
tion with its amplitude decreasing. Hence, the dynamics of the boost-invariant 
magnetic field is stable unless one allows for rapidity-dependent fluctuations. 
The dissipative behavior of the solution is understood as the effects of longitu- 
dinal expansion of the system. If the second term in Eq. (23) were absent, the 
solution would be simply given by a non-dissipative oscillation e*' fc± ' T . Thus 
the dissipative behavior is a result of the second term, which is present due to 
the r-rj coordinates that are suitable for describing expanding systems. 

To get an intuitive picture of the time dependence of the solution (24), let 
us discuss a simplest but physically reasonable setting for the initial condi- 
tion. Consider a single isolated magnetic flux tube which is represented by a 
Gaussiarf 3 ! with its transverse size being 1/Q S - This can be realized by the 
following initial condition: 

S^JfcjJcx-z^e'Sr (26) 
k ± 

In fact, one can compute B z for r > with the result 

B*( T , r) = B e- Q " r2 e- Q " T2 I {2Q 2 s rr) , (27) 



where r = \x±\ and Iq(z) is the modified Bessel function. In this expression 
one can easily check that B z (t = 0,r) is the Gaussian. Similarly, one can 

compute the transverse electric field E T = ^ (E x ) 2 + (E y ) 2 as 

E t (t, r) = B e -^V Q * r2 h(2Q 2 s rr) . (28) 

3 We can perform similar analyses for Lorentzian profile l/{k \ + Q 2 ) in the mo- 
mentum space. 
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Fig. 2. Spatial profile of the longitudinal electric field B z (left) and the transverse 
magnetic field E T = a/ (E x ) 2 + (Ey) 2 (right) at five different times Q s t = (solid), 
0.5 (longest dash), 1.0, 1.5, 2.0 (shortest dash). The maximum strength of the 
electric field is normalized to unity at r = 0, i.e., B = 1. Notice that E T = at 
r = 0. These are true for E z (left) and B T (right) of a single electric flux tube. 

In Fig. 2, we show the transverse profile of B z and E T at five different instants 
Q s t = 0, 0.5, 1.0, 1.5, 2.0. Remarkably, as r increases, the magnetic flux tube 
expands outwards while the strength inside the flux tube decays rapidly. On 
the other hand, the transverse component of the electric field is initially zero, 
but grows until Q s r ~ 1 and then turns to slow decrease. 

All these features are intuitively understood by looking at the evolution of 
the flux tube as a function of t = x°, instead of r. In Fig. 3 shown are the 
lines of the magnetic field in the z-x± plane at two different instants Q s t = 1 
(solid) and 2 (dashed). The field strengths at Q s z = ±1 (±2) for Q s t = 1 
(2) have the Gaussian profile of the initial condition. Non-zero transverse 
magnetic field appears at z ^ 0, according to Eq. (19). As the nuclei recede 
from each other, the flux tube in between is longitudinally stretched out and 
transversely expands, and thence the field strength decays rapidly in the mid- 
rapidity region. Notice that this behavior is completely the same as that of 
ordinary electric field since we have assumed that the gauge field has a single 
color component. 

It is interesting to note that the rapid decay of the longitudinal magnetic field 
and the slow growth of the transverse electric field is qualitatively similar to 
the numerical result reported in Ref. [6]. In order to confirm this similarity in a 
more direct way, let us take the average of (B z ) 2 and (E ) 2 over the transverse 
plane. For a single isolated flux tube given by Eq. (39), the integration over 
the transverse plane gives finite results: 




2 Ql 



UQlr 2 ) , 



(29) 
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Fig. 3. Lines of the magnetic force at Q s t = 1 (solid) and 2 (dashed). The sources 
are located at the longitudinal positions, Q s z = ±1 and ±2, at Q s t = 1 and 2, 
respectively. These can be also regarded as the lines of the electric force in the 
electric flux tube. 




Fig. 4. Time dependence of the magnetic field (B z ) 2 (solid) and the transverse 
electric field (E T ) 2 (dashed) of a single magnetic flux tube. Both are averaged over 
the transverse plane, and normalized for (B z ) 2 at r = 0. These are true for (E z ) 2 
(solid) and (B T ) 2 (dashed) of a single electric flux tube. 

W? = j d 2 x ± (E T ) 2 = \ *Le-^h(Qy) . (30) 



In Fig. 4, we show the time dependence of (B z ) 2 and (E T ) 2 which is nor- 
malized for the magnetic field at r = 0. This is qualitatively very similar to 
the numerical result of Ref. [60- For more realistic initial configurations with 
many flux tubes (of electric and magnetic types), we have to sum up all those 
contributions. However, if the flux tubes are distributed homogeneously on 



There is no divergence in our computation since we use a Gaussian profile in the 
transverse plane for the flux tube [14] . 
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the transverse plane, which will be realized for the initial condition of well- 
saturated nuclei, we can expect that the contribution from a single flux tube 
gives a fair approximation to the actual behavior of (B z ) 2 and (E T ) 2 aver- 
aged over the transverse plane. Hence it makes sense (at least qualitatively) 
to compare our results with the numerical results of Ref. [6]. 

The large-r behavior of (B z ) 2 and {E T ) 2 can be easily determined from the 
asymptotic form of the modified Bessel function: 

e z 

I v (z) ~ — , z — > oo . (31) 



Namely, both (B z ) 2 and (E T ) 2 show the same behavior as is evident from the 
figure: 

JWf ~ JWf ~ - . (32) 

T 



This is consistent with the scaling behavior of the energy density at large r, 
which is a result of the longitudinal expansion of the system [14]. Indeed, as 
we will show later for the electric flux tubes, the behavior (32) is common for 
the longitudinal electric field and transverse magnetic field. 



These time dependences of (B z ) 2 and (E T ) 2 are also understood from the 
viewpoint of the energy conservation. The problem can be treated just as 
in the case of the 1-dimensional Bjorken expansion since we have integrated 
the density over the transverse plane. Change of the energy density per unit 
rapidity is caused by the pdV work [21], where p is the longitudinal pressure 
and dV is the volume change dV = dr per unit rapidity for a boost-invariant 
expanding system. In the magnetic flux tube configuration, the energy density 
and the longitudinal pressure [7] integrated over x± are, respectively, given by 
\[{E T ) 2 + {B Z ) 2 \ and \[(E T ) 2 - (B z ) 2 ] at rj = 0. Then the energy conservation 
implies 



d 



' T 



(e t ) 2 + (B z y 



{E T ) 2 - (B 



z\2 



dr 



(33) 



As is well known, without the pressure term, the energy density per unit ra- 
pidity would decrease as 1/r. However, for the initial configuration having the 
longitudinal field alone, the right-hand side is positive [7,6], which slows down 
the decrease of the energy density. As the transverse field (E T ) 2 is generated, 
the longitudinal pressure is weakened and the behavior of the energy density 
approaches the 1/r dependence. 
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2. 3 Electric flux tube 



Let us now consider the opposite situation where only a longitudinal electric 
field is present. This situation is simply realized, from Eqs. (16) and (17), 
when a\ and a 3 2 (i ^ j) have the same color but a\ and a\ have different 
colors from each other. In this case, «j is pure-gauge in the transverse plane 
and can be eliminated by a gauge transformation. Therefore we consider the 
configuration 

a ± 0, OLi = 0. (34) 



Then, each component of the field strength is given by (at r\ = 0) 



E? = 0, E z = -d T (—T 2 a) , (35) 

T 

B i = reijdja , B z = . (36) 

We again assume that a is oriented to a certain color direction. This color 
direction of the isolated electric flux tube is independent of that of a 1 for the 
magnetic flux tube in the previous subsection. While the second and third 
Yang-Mills equations, (12) and (13), are automatically satisfied under this as- 
sumption, the first equation (11) reduces to a very simple differential equation: 



-d T (r 3 d T a) -dfa = 0. (37) 



Performing the Fourier transform in the transverse space, one finds that this 
equation is identical with the Bessel equation. Therefore, the solution satisfy- 
ing the initial condition a(r = 0,x±) = ai n i t (x±) is given by 

a{T,x±) = I ^ elfcxx± Fk Jl(|A;±|r)5init(A;±) ' (38) 



where (5mit(^±) is the Fourier transform of the initial condition ai n a(x±_). 
One can easily check that this solution satisfies the initial condition by using 
lim T ^ jj^Ji(\kj_\T) = 1. Similarly as in the case of the magnetic flux tubes, 
the expansion of Ji ( | A:_i_ | t) / 1 | r with respect to r gives only even powers of 
r, which is consistent with the result of r-expansion developed in Ref. [20]. 
In addition to this, we note that J\(z)jz is an oscillating function with its 
amplitude decreasing. Thus, the dynamics of the boost-invariant background 
electric field is stable, too. 
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Let us again consider a simple example in order to understand the time depen- 
dence of the solution (38). Consider a single isolated electric flux tube whose 
transverse size is of the order of 1/Q S . This can be represented by the initial 
condition of the Gaussian form: 

a- m \t(x±) = a e~^ Q * . (39) 

Inserting (the Fourier transform of) this into Eq. (38), one finds that the result 
is represented by infinite series of hypergeometric functions. However, the 
expressions for the longitudinal electric field E z and the transverse magnetic 
field B T = y/ (B x ) 2 + (By) 2 become remarkably simple: 

E z {t, r) = E Q e- Q > 2 e~ Q ' T2 1 {2Q 2 s rr) , (40) 
B t (t, r) = E e -Q^\-Q 2 r 2 l^2Q 2 a rr) . (41) 

Notice that these are the same functional form as B z and E T in the previous 
subsection. Thus we can simply replace all the results discussed for B z and 
E T by E z and B T . In particular, Figures 2 and 4 are equally true for E z and 
B T , and Figure 3 is true for the electric field. 

2.4 Coexisting case 

So far, we have discussed magnetic and electric flux tubes separately. However, 
in fact, if both a and a, are oriented to the same direction in the color space, 
then one can ignore all the commutators in the Yang-Mills equations (11)— (13) 
yielding decoupled equations for on and a. Namely, they are given by Eq. (21) 
and Eq. (37), respectively. Moreover, as is evident from the explicit form of 
the field strength (5)-(8), the contributions to the field strength from a and 
a, are additive. Namely, for 77 = 0, one obtains E l = d T at, E z = ^d T (—r 2 a), 
B % = reijdja, and B z = —d x a y +d y a x . Therefore, all the results in the previous 
subsections are valid even in the case where both a and are nonzero as far 
as they have the same color structure. 



3 Instabilities of the glasma 

In the previous section, we have studied the dynamics of boost-invariant 
glasma as a collection of flux tubes which have typical diameters of 1/Q S and 
extend between the two receding nuclei. There we have learned that the flux 
tubes (longitudinal electric and magnetic fields) expand towards outside in 
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the transverse directions, and the strength inside the tubes rapidly decreases. 
On the other hand, the transverse components of the field strength B T and 
E T are zero at r = but grow slowly until around Q s t ~ 1, where the mag- 
nitudes of B T and E T become comparable with the longitudinal components 
E z and B z . This fact, however, does not imply isotropization of the system. 
The decay of the longitudinal fields is, in our picture, due to the transverse 
expansion of the flux tube, which also induces nonzero transverse components 
of the field strength. Thus, all these can be understood as the process towards 
homogeneity in the transverse plane. In fact, the gauge fields a and CKj are 
still rapidity independent. Therefore, thermalization or even isotropization of 
the system never occurs with this boost-invariant glasma. This is where the 
rapidity-dependent small fluctuations come into play. 

In this section, we discuss the stability of boost-invariant glasma against 
rapidity-dependent perturbations. The small fluctuations are defined by 

Ai = Ai + ai(r, rj, x±) , A v = A v + a v (r,r],x±) , (42) 



where the background fields Ai and A v are boost-invariant glasma discussed in 
the previous section (We stay in the Fock-Schwinger gauge A T = A T + a T = 0). 
The stability of the boost-invariant glasma can be examined by analyzing the 
Yang-Mills equations (l)-(3) to the linear order in the fluctuations a* and a n 
on the background fields. The non-Abelian nature uniquely fixes the coupling 
between the background field and the fluctuations. 

In our calculation, to make the situation simple, we approximate the boost- 
invariant background fields by transversely constant electric or magnetic fields 
which are of course rapidity independent. This simplification is a reasonable 
approximation for the fields inside the flux tube. Thus, we implicitly assume 
that the transverse extent of the constant electric or magnetic field is of the 
order of l/Q s (or slightly larger than that), and we will later examine if the 
typical transverse size of the unstable mode is smaller than l/Q s . In studying 
the equations for the fluctuations, we restricted ourselves to the £77(2) color 
gauge group and greatly benefited from Ref. [22] for technical details. 

3.1 Fluctuation in the magnetic background field 
3.1.1 Constant magnetic field 

Consider the simplest case where only time-independent and spatially constant 
magnetic field is present as the background field. This is a special case of 
Sect. 2.2 where we treated only CKj(r, x±). In the following, we consider color 
SU (2) group, and assume that the constant magnetic field is directed to the 
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third color component. This is realized by 

a = 0, a? = 5 a3 ^(y5 tl -x5 t2 ), (43) 

where B(> 0) is independent of r. The equations for fluctuations are coupled 
equations for and Oj, and are still complicated as they are. If one assumes 
that a v 7^ and a« = 0, then one can show that a v is stable. We expect this is 
true even in the full coupled equations, and we rather consider the opposite 
case: a T] = and a, ^ . Then the Yang-Mills equations are 



[Di, d v ai\ = , (44) 

[T>i, rd Tai ] = , (45) 

— d T {Td T ai) + — dldi + [Vj, SFji] - %g\a^Tj^ = , (46) 

T T z 

where SJ-'ji = [Dj, Oj] — [X>j,Oj]. The 77 dependence can be easily treated by 
introducing the Fourier transform. The first equation implies: 

[2^1 = 0. (47) 



This is the Gauss law. If this is satisfied, the second equation is trivially satis- 
fied, too. The Gauss law is still very helpful in simplifying the third equation. 
After a long manipulation, we can rewrite the third equation as follows: 

--d T (rdX) ~ ^~tf + dial + gBe Zah d B ~a\ 

T T 

+2gBe 3ab (a b 2 5 tl - ~a\5 i2 ) + ^^ e 3cb e 3ba ~a<r = , (48) 
where we have introduced the Fourier transform with respect to r\ 

o*(t, 77, x±) = J ^ e™ v di(r, v, x±) , (49) 

and cylindrical coordinates 

x = rcos#, y = rsin9. (50) 



At this stage, it is obvious that the 'neutral' fluctuation a 3 is stable. The 
solution of Eq. (48) without the coupling to the constant B, is given by the 
Bessel function a 3 oc Ji U (\k±\r), which shows a dissipative oscillation. 
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On the other hand, the other color components are not diagonal, and couples 
with each other through the magnetic field B. Physically, these fields behave 
as 'charged' matter coupled to the gauge field, and show cyclotron motion. 
Thus it is convenient to introduce eigen-states of angular momentum operator 
in the transverse plane, which is given by the ^-derivative: 

d 

— = xd 2 - ydi = -iL . (51) 



Taking the eigen-states of L, i.e., af oc e ime , we can replace dg by im. Then 
we find that the transverse structure is equivalent to the two-dimensional 
harmonic oscillator. This is clearly seen if we introduce the combined fields 
for colors 

af } = ~a\±i~a\, (52) 

and moreover for the transverse coordinate 

a£ = a[ } ± i a { 2 } . (53) 



Now the equation is diagonal with respect to these fields: 



\d T {Td T ~a¥) + (^ T mgB ± 2 9 b} af> + (-d\ + a<? = , 

l -d T {rdJ ±] ) +{^ 2 T mgB T 2gB^ a™ + + 9 -^^j = . 

The last terms can be identified with the Hamiltonian of the two-dimensional 
harmonic oscillator. Therefore, one can replace them by the energy eigen- 
values E n a^ with E n = (2n+\m\ + l)g B . Then we finally obtain the equations 
of the familiar form: 



-d T (Td T af>) + \^E n T mgB ± 2gB + - af> = , (54) 

\{rd T ~a { ^) + {E n T mgB t 2gB + ^ a {±) = . (55) 

These are equations for the Bessel functions. Notice that the sign of E n =f 
mgB ± 2gB or E n =p mgB =p 2gB crucially affects the behavior of the solution. 
If it is negative, the solution is unstable. For the + component of the spatial 
coordinates (Eq. (54)), it can be negative for n — 0, color — , and m < 0, 
i.e., Eq + mgB — 2gB = —gB. On the other hand, for the — component of 
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the spatial coordinates (Eq. (55)), it can be negative for n — 0, color +, and 
m > 0, i.e., E —mgB — 2g B = —gB. Both cases reduce to the same equation 
which has an unstable solution: 

dlf+ l -d T f+(-gB + y\f = Q. (56) 



The solution is given by the modified Bessel function: 

a$r\ aL +) oce' m VW e -^/ w (^ T ) . (57) 



These are the exact solutions to the linearized equations for the fluctuation 
Oi(r, 7), x±) in a constant magnetic field. The functional dependence on m 
and r = comes from the wavefunction of the two-dimensional harmonic 
oscillator with n — 0. The r dependence is completely determined by the 
modified Bessel function, which is consistent with the results suggested in 
Ref. [14]. By using the asymptotic form (31) of I v {z) whose leading term 
is independent of v, we indeed find that the above solution is exponentially 
divergent: 



s(-) ~(+) e 



a>->, a>y ~ , = (58) 
2-k^/JBt 



Therefore, instability exists for any v although appearance of the instability 
depends on v, as we will discuss later. All the other modes except for those 
mentioned above are stable, and their time dependence is governed by the 
Bessel function J iv {^fgBT). 

As we already mentioned, the fluctuations as the charged matter obey 
cyclotron motion in a uniform magnetic field. The unstable modes appear in 
the lowest Landau level (n = 0). Interestingly, the fluctuations o^', having the 
opposite sign of charge, show instabilities in the opposite cyclotron rotations 
to each other. In fact, the existence of instabilities in a constant non-Abelian 
magnetic field has been known in the ordinary cartesian coordinates [22,23]. 
What we have done here is essentially to find the same instability phenomena 
in the r-rj coordinates (see Ref. [16] for a similar conclusion). 

Let us investigate more about the properties of unstable modes. First of all, 
the smallest spatial size of the instability is given by m = 0. Namely, the 
Gaussian profile in the two dimensional plane has the spatial size L given by 
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®e{I iv (z)} 




On the other hand, the unstable mode with \m\ ^ is distributed in a ring of 
radius £ ~ ^2\m\/gB. For the magnetic field created in the CGC, gB ~ Q 2 S , 
one finds L ~ 1/Q S . Recall that the transverse size of the flux tube is ~ l/Q s 
at t = 0, but it will expand outward in the transverse plane. Therefore, we 
conclude that the magnetic instability (57) will indeed be realized in the time 
evolution of the glasma. 

Next, we carefully investigate the time evolution of the unstable modes (57) 
which is completely determined by the modified Bessel function I{ V {z). In 
Fig. 5, we show typical behaviors of $te{Ii„(z)} for two different values of v. 
The figure indeed shows divergent behavior suggested by the asymptotic form 
(Eq.(31)). The time scale for the instability to grow is given by 



1 



1 



' grow 



(60) 



However, the time when the instability sets in depends on the value of u, which 
we denote as z wait (because we have to wait until z = z wait ). From Fig. 5, one 
may infer a rough estimate 

^wait ~ V ■ (61) 



This can be indeed justified in two different ways as we discuss below. Consider 
first the explicit form of Ii V (z): 
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Notice that the term (z/2) %v = e lulnz / 2 simply oscillates for real is and z, and 
thus it is not relevant for the instability. One can estimate z wa i t from the other 
part. Taking the first two terms in the infinite series, one finds 



T(iv + 1) r(iz/ + 2)J T(iz/ + 1)[ 1 + z/ 2 V2 

Therefore, one can naively define the position z WSuit where the second term 
becomes the same order as the first one. For the real part of Ii v (z), this yields 



V^wait/2) 



2 



k (63) 



for k 2 = 0(1). Therefore, for large is, z wa n grows linearly as a function of is: 



"wait 



2«Vl + is 2 ~ 2ku. (64) 



Choosing 2k ~ 1 leads to the rough estimate (61) inferred from Fig. 5. 

Another way to estimate z wa it is based on the differential equation. The mod- 
ified Bessel function I iu ((3z) is a solution to /" + f'/z + (— (3 2 + is 2 /z 2 )f = 0. 
The behavior of the solution change depending on the sign of the coefficient 
of the last term — (3 2 + v 2 / 1 z 2 . If it is positive, which is realized at small z, the 
solution is oscillating. On the other hand, if negative, then unstable. The sign 
change takes place at — [3 2 + is 2 / z 2 = 0, yielding the same estimate for z wait , 
i.e., z wait = v for (3 = 1. 

If we apply this notion to our result, we can estimate the time r wa ; t when the 
instability sets in (which may be called "waiting time"). For z = \J gB r ~ 

Q s r, 

QsTwait ~ is . (65) 



Therefore, for a fixed value of is, the instability will exist only for r ^ T wait . 
In fact, this waiting time Q s r wa i t can be much longer than the growth time 
Qs^grow ~ 1 for large is. Although the instability grows with the characteristic 
time given by r grow once it starts, we have to wait until r wait . Therefore, the 
total time for the instability to become manifest is estimated as 

<2 S (r W ait + r^w) ~ V + 1 • (66) 



It is essentially determined by the waiting time r wait for large is. But this time 
scale can be short for small values of is. Since the plasma instability scenario 

5 We thank A. Iwazaki for this point. Similar argument can be found in Ref. [16]. 
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is based on the kinetic approach which is valid only for Q s r ?t 1, we expect 
that the glasma instability will be more operative for the rapid thermalization 
of the matter created in the heavy-ion collisions. 

The fact that the mode of v has the "waiting time" r wait for the instability 
implies that at fixed r there exists a maximum value of v which can participate 
in the instability, and this maximum value increases with increasing r. In the 
present case, it is given by 

^max ~ QsT • (67) 



The existence of z/ max and its linear dependence on r are consistent with the 
observation numerically found in Ref. [7]. But the relation between z/ max and 
r, such as the power of r or the coefficient, crucially depends on the behavior 
of the background field as we will discuss below. 



3.1.2 Time- dependent magnetic field 

So far, the longitudinal magnetic field was assumed to be time-independent, 
but the basic equation (48) is in fact true even for time-dependent background 
field. However, we have to be careful not to replace simply the harmonic 
oscillator part by its eigenvalues when the magnetic field is time-dependent. 
This is physically because the cyclotron radius changes in accord with B(t), 
which invalidates the factorization of the transverse dynamics from the time 
dependence. Still, we may assume a situation that the time-dependence of 
the magnetic field is slow compared with the characteristic time scale of the 
fluctuation (namely, the growth time of instability), which corresponds to 
neglecting time derivative of B{r). Thus we can recover the same equation as 
before except that the magnetic field B is replaced by the time-dependent one 
B(t). For example, the lowest Landau level yields the following equation 

d 2 T f + \d T f + (-gB(r) + ^) f = . (68) 



This equation can be solved for the following simple case: 

B(r) = B . (69) 



The solution is again given by the modified Bessel function (for 7 = 2, the 
solution does not show rapid growth. It is given by a power of r, or just 
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oscillating): 

/M-^^^r) , P=^, 7^2. (70) 

As we mentioned before, the relation between z/ max and r depends on the 
precise form of the background magnetic field: 



iW ~ yJgB(r)T. (71) 

Linear dependence appears only for a time-independent magnetic field as seen 
in Eq. (67). 

For the special case of 7 = 1, one finds 

/ ~ I 2w (2yfgBoT S) (72) 

which is similar to what is expected for an expanding system. On the other 
hand, if one takes 7 = 1/2 which is consistent with the asymptotic behavior 
(32), the solution is given by 

/ ~ W (y 9B To /2 r 3/4 ) , (73) 
which is consistent with the result of Ref. [16]. 

For a generic time-dependent magnetic field B(t), the equation is no longer 
solved by simple i^'s but must be solved numerically. Still one can infer the 
qualitative behavior of the solution from the third term of Eq. (68). Assume 
that B(t) > is slowly decreasing from a certain finite initial value B{r = 0). 
Then, the coefficient of the third term (—gB(r) + v 2 /t 2 ) is positive for small 
r, but will turn to be negative at a later time. Therefore, the solution is 
oscillating until r = r wait determined by — g-B(r wai t) + ^ 2 / r wait = 0, and then 
starts to grow, which is qualitatively in agreement with the behavior of the 
modified Bessel function shown in Fig. 5. 

Notice that the result (71) for z/ max is still valid. In particular, if one substitutes 
Eq. (29) as the time-dependent magnetic field, one obtains (B = B /2) 

1/4 



gB (e~^ 2 1 (Q 2 r 2 )) ' r. (74) 



As the factor (e °= r2 /o( < 5s r2 )) ^ decreases quite slowly, one may observe an 
approximate linear increase of z/ max in r, in accordance with the numerical 
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result [7]. But more precisely, the factor depends on r as ~ 1/r 1 ^, and thus 
^max ~ t 3 ^ 4 as already suggested in Eq. (73). 

3.2 Fluctuation in the electric background field 
3.2.1 Constant electric field 

Consider the case where only time-independent and spatially constant longi- 
tudinal electric field is present as the background field. This is a special case 
of Sect. 2.3 where we treated only a(r,x±). Indeed, a constant electric field 
E z = E is realized by the following gauge field: 



Although the equations for the fluctuations a v and a, in this background are 
coupled equations, we can simplify them similarly as in the magnetic case. If 
one assumes that a T] ^ and a, = 0, then one can show that a v is stable. 
We expect this is true even in the full coupled equations, and we consider the 
opposite case: a v = and ^ . The Yang-Mills equations are now simplified 
to the followings: 



In order to solve these equations, it is convenient to introduce the Fourier 
transform of with respect to i] and x± 



a 




a* = 0. 



(75) 




(76) 
(77) 

(78) 




(79) 



and its perpendicular component b via 




(80) 



Then, the third equation can be rewritten as 



-d T (Td T ~b) + k 2 ± b + 




(81) 
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while the first and second Yang-Mills equations are automatically satisfied by 
the introduction of b (because fcjOj = 0). 

In order to further simplify the equation, let us consider the SU (2) case and 
assume that the electric field is oriented to the third color direction, i.e., 
E = E a T a with E 1 = E 2 = and E 3 = E. Then, one obtains 

1 1 / ~ n 2 E 2 T 4 \ 

-d T (rd T b a ) + — v 2 b a - iugET 2 e 3ba b b + ^ e 3ca e 3cb b b + k]b a = . 

r T \ 4 J 

For the 'neutral' fluctuation a = 3, the equation does not depend on E, and 
the solution is given by the Bessel function b 3 oc Ji U (\k±\r). Therefore, there 
is no instability in the fluctuation to the third color direction. On the other 
hand, the equation is not diagonal for a = 1,2. One can easily diagonalize it 
by introducing 

~b {±) =~b l ±ib 2 . (82) 
Then one finds 



(d 2 + \) 6<±> + jjfci + 1 [v 2 ± ugEr 2 + ^r 4 ^) | 6<±> = . 



(83) 



Let us first ignore the third term in the bracket since it is proportional to r 4 
and is small enough compared to the other terms at the earliest times. This 
approximation should be valid for r satisfying 



For the electric field in the CGC gE ~ Qj, this condition reads 

Q s r<^V2^. (85) 

Within this approximation, one obtains a simple differential equation: 

d 2 T + -d T ) 6 {±) + ( k 2 ± ± ugE + 41 & (±) = . (86) 



t J \ r 2 



This is nothing but the equation for the Bessel functions. As we discussed in 
the magnetic case, the behavior of the solution depends on the sign of k\±vgE. 
If we assume that both E and v are positive, then, for the plus sign, it is always 
positive, k\ + vgE > 0, and the solution is given by oc J; w {\jk\ + vgE r). 
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Therefore, the fluctuation is stable. On the other hand, for the minus sign, 
if k\ — vgE < 0, then the solution is given by the modified Bessel function 



U-\t, u, k ± ) oc IMugE -klr). (87) 



As has already been seen, the modified Bessel function Ii V (z) shows unstable 
behavior for large value of z, which corresponds to late time. We shall check 
below whether such an instability will indeed manifest itself within the time 
specified by Eq. (85). 

The condition k\ — vgE < can be regarded as the condition for the existence 
of unstable fluctuations. In other words, unstable fluctuation has transverse 
extent L ~ l/k± larger than 1 / -JvgE. Recall that the transverse size of the 
coherent gauge field is of the order of l/Q s . Therefore, unstable fluctuations 
may exist if the transverse extent L is smaller than l/Q s : 



Q s \JvgE y/uQ s ' 



This is satisfied if the momentum v conjugate to rapidity is large: 

v > 1. (89) 



Since the solution is again given by the modified Bessel function Ii U (z), all 
the properties discussed for the magnetic instability can be applicable to the 
present case. Notice that both the waiting time r wa i t and the growth time r grow 
for the instability depend on k"j_, and that the shortest times are realized for 
k± = 0. Namely, for fixed u, the shortest waiting time and the shortest growth 
time are, respectively, given by 



Qs^wait^V^, (90) 

<2 s r grow ~ \/y/v. (91) 

If the waiting time Qs^wait plus the growth time Q s T grovi is within the validity 
region (85), one may be able to conclude the existence of instability. However, 
this condition reads 

Q s {r waii + T- grow ) ~ ( 1 + -) < V2^ , (92) 



which is, unfortunately, not satisfied by any value of v > 0. Therefore, the 
approximate equation (86) allows for unstable solution, but it shows the in- 
stability only outside the region of validity. 
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Indeed, this means that one cannot ignore the the E 2 term in Eq. (83), and 
then one finds that its solution does not show instability. This is immediately 
understood if one notices that the second term in the curly bracket can be 
written as [y ± gEr 2 /2) 2 /r 2 which is always non-negative. Thus, the curly 
bracket in Eq. (83) never becomes negative. On the other hand, the instability 
in the approximate equation (86) occurs when the corresponding part becomes 
negative, which is not allowed in the original equation (83). And the upper 
limit of the time (85) corresponds to the time when [y — gEr 2 /2) 2 jr 2 = 0. 

In fact, one can write down the explicit form of the solution to the original 
equation: 

oc \ [c + M K „ (ifr 2 ) + c_M K ,_, (ifr 2 )) , (93) 



where M Kitl (z) is the Whittaker function and k and /i are both pure imaginary 
numbers 

i ( k 2 , \ v , , 



By using the asymptotic form of the Whittaker function, one finds that there 
is no instability at large r. 

What we have done should be essentially equivalent to the Schwinger mecha- 
nism. Unlike the creation of massive fermion and anti-fermion pairs, there is no 
mass gap for the 'charged' gluon matter &W. Thus any external electric field 
allows pair creations. In fact, the 'charged' fluctuations are accelerated in the 
electric field, and acquire large longitudinal momentum u±gEr 2 /2 = u±gA Tj 
[24]. Therefore, even if the fluctuation does not grow in time, it will be in- 
finitely accelerated towards the ends which are moving at the speed of light. 
This situation may be also called instability as discussed in Ref. [22]. It should 
be possible to find more explicit correspondence between our case and the 
Schwinger mechanism, and it will be interesting to also consider the effect of 
backreaction. We leave such topics for future study. 



3.2.2 Time- dependent electric field 

Let us briefly discuss the case with time-dependent electric field. Although 
we can in principle treat any time- dependent electric background, we here 
discuss only the following simple case where we can easily find the solution to 
the equation: 

E(r) = E (^y, 7 ^2. (95) 
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To obtain the differential equation for the fluctuations, one has to replace E 
in Eqs. (83) and (86) by 2£(r)/(2 - 7) = (5E{t) (instead of E(t)). Then, for 
the approximate equation (86), 




6^ = 0. 



(96) 



One can reproduce the time-independent case by setting 7 = (a = 1). 
When k± = 0, the solution is again given by the modified Bessel function. For 
example, 



This solution should be a good approximation for small r, but the unstable 
behavior at large r is just an artifact which appears when we ignore the E 2 
term. 



4 Summary and discussion 

We have investigated the dynamics and stability of the glasma based on the 
picture that the glasma is initially made of longitudinally extending flux tubes 
with their transverse size being 1/Q S . For simple Gaussian flux tubes where 
the longitudinal fields have only one color component, we found that the flux 
tubes expand outwards as time goes, and the transverse field strengths are 
accordingly induced. This time evolution is quite consistent with the numer- 
ical result reported in Ref. [6], and thus we expect that our simple picture 
captures the essence of the actual physics situation. It is very remarkable that 
the Abelian dynamics of the background field in our model reproduces the 
qualitative behavior of the time-dependence of the energy density found in 
Ref. [6]. At the same time, this raises a question what is the role of the non- 
Abelian nature of QCD in the result of Ref. [6] . One of the possible effects of 
non-Abelian dynamics will be the interactions between flux tubes with differ- 
ent colors. This was ignored in our calculation for an isolated flux tube, but we 
expect such effects will become relevant only at later time when the overlap 
between the expanding flux tubes is sizable. 

We have studied in detail the response of small rapidity-dependent fluctua- 
tions in the background of boost-invariant glasma. Instabilities in the electric 
or magnetic fields are, if any, expected for the 'charged' fluctuations since 
they can couple to the background field. Recall that there is no preferred 
color direction in the Yang-Mills theory. Thus, if the field is forced to have 




(97) 
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only a certain color component, it is quite natural that fluctuations perpen- 
dicular to the forced color component (namely, the 'charged' fluctuations) 
become large, and can grow exponentially in some case. In fact, we found 
that the 'charged' fluctuations can grow exponentially only in the magnetic 
background. On the other hand, the response of the 'charged' fluctuations in 
the electric background is essentially equivalent to the Schwinger mechanism. 
While the 'charged' fluctuation is accelerated infinitely and the system is un- 
stable in the Schwinger mechanism, the fluctuation itself does not grow in 
time. 



The unstable mode in the magnetic background corresponds to the lowest 
Landau level whose transverse size is 1/Q S . Thus it can survive in the flux tube 
which expands in the transverse directions with its initial size being 1/Q S . The 
time for the instability to become manifest is estimated as Q s r ~ v + 1 where 
v is the momentum conjugate to the rapidity coordinate. Thus, for small v, 
instability may occur at very early stage of the collisions, even before the 
ordinary plasma instability starts to operate. Thus, we hope that the glasma 
instability play an important role in the early thermalization. 

As we have already emphasized, the coupling between the 'charged' fluctuation 
and the magnetic background field is the essential dynamics for the existence of 
instability. Since the 'charged' fluctuation itself is the gauge field, such coupling 
is possible only in the non-Abelian Yang-Mills theories. Therefore, one may 
say that the glasma instability is a result of non-linearity of the Yang-Mills 
equation. Having said this, we recall the similarity of the glasma instability 
with the emergence of chaos in the classical Yang-Mills theory [25] . It is known 
that the spatially constant Yang-Mills fields exhibit chaos due to the non- 
linear interactions among the fields. Notice that if there are unstable modes 
as in the glasma, the Lyapunov exponent for two unstable modes with slightly 
different initial conditions is given by the growth constant of the unstable 
modes. This suggests an interesting viewpoint of the Yang-Mills chaos to the 
glasma instability. 

The Weibel instability has been widely studied as the main mechanism of 
the plasma instability in the context of heavy ion physics [3]. Although it is 
not evident at present whether the magnetic instability found in the present 
paper is equivalent to it, our setup is quite similar to that of the plasma 
instability scenario. For example, the 'hard' particles correspond to the boost- 
invariant background field in our study, which are supplied by the CGC initial 
condition, and the 'soft' fields are the rapidity-dependent fluctuations in our 
study. Since the background field is independent of rapidity, it yields a very 
anisotropic distribution in the momentum space (k± ~ Q s , k z ~ 0). Moreover, 
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linearizing the Yang-Mills equations for the fluctuations a 



K^[A]a u = 0, 



is nothing but the computation of the inverse propagator K^ U [A] of the 'soft' 
fields on the background field A^. This may correspond to the hard- loop cal- 
culation in the study of the plasma instability. We expect that future investi- 
gation will clarify the possible relationship between our result and the Weibel 
instability. 
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